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(54) Bandwith estimation and adaptive filtering 

(57) The present invention concerns a method for 
bandwidth estimation and a method for adaptive filtering 
in a data receiver, wherein said data receiver has been 
provided with a set of filters each having a known band- 
width and capable of removing noise from a quantity z r 
derived from a received signal, either method compris- 
ing, for each one of these filters: a) the application of the 
filter to z r , thus yielding a quantity z F , and b) the calcu- 
lation of the power A of the difference between z r and 



z F . For bandwidth estimation, one determines for each 
filter whether it has truncated the noiseless component 
of z r , on the basis of a criterion involving A and a 2 n , 
where a 2 n is an estimate of the power of the noise com- 
ponent of z r . For adaptive filtering, one selects one filter 
out of the set on the basis of the same criterion, in such 
a way that the largest possible amount of noise is re- 
moved from z r without truncating its noiseless compo- 
nent. 
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Description 

[0001] The present invention concerns the presence of additive noise in received signals. More specifically, it con- 
cerns the estimation of the bandwidth of the noiseless component of any chosen quantity derived from a received 
5 signal (possibly this received signal itself). It also concerns the adaptive suppression of the noise component of such 
a chosen quantity. 

[0002] The invention is particularly useful in the field of wireless digital communications systems. In these systems, 
one usually deals with multipath channels: due to reflections on obstacles, the transmitted signal reaches the receiver 
after having followed a plurality of paths. Hence the received signal results from the superposition of a plurality of 

10 replica of the transmitted signal, each one associated (in the equivalent complex baseband representation) with a 
specific delay and attenuation. This received signal is therefore equal to the convolution of the transmitted signal with 
a so-called Channel Impulse Response (CIR) h(f) (where t denotes time). Equivalently, implementing a Fourier trans- 
form, this time-dispersive characteristic of a channel may be represented in the "frequency domain" by means of a 
"channel transfer function" H(f); at any given frequency, the value of this function is called the "Channel Transfer Factor" 

15 (CTF). The CIR is non-zero only in an interval 0 < t < t, where t is the so-called excess delay of the channel, viz., the 
maximum possible delay between the various paths. 

[0003] The data symbols to be transmitted are modulated, and the receiver demodulates the corresponding received 
signal in order to recover those symbols. Clearly, knowledge of the CIR allows one to conduct this demodulation more 
accurately. 

20 [0004] Furthermore, in many modern digital communications systems, for example telephone systems or mobile 
radio systems, the demodulation process is "coherent". This requires that the system provides a means for "channel 
estimation", viz., for computing an estimate of some channel characteristics of interest (such as the channel impulse 
response) using a reference portion of the received signal. For example, in so-called pilot-assisted channel estimation, 
one transmits a "pilot signal" based on symbols which are a priori known by the receiver. As another example, in so- 

25 called decision-directed channel estimation, the receiver uses as a reference some received symbols whose value 
has been determined by means of a tentative decision process. 

[0005] But a realistic channel always suffers from random noise, which adds up to the ideally-noiseless received 
signal. Therefore, if channel estimation is purely based on reference symbols, one only obtains a coarse channel 
estimate, so that the coherent demodulation process based on this estimate will offer sub-optimal performance. In 
30 order to improve over such coarse estimates, it is often necessary to take into account known statistical properties of 
the channel and/or the noise process. For example, channel variations over time or over frequency have typically a 
smaller bandwidth than the noise process, so that the amount of noise afflicting the estimates can be reduced by using 
a filter. 

[0006] A filter is a device for signal processing whose output, in the frequency domain, is equal to the product of its 
35 input and of a complex frequency-dependent gain factor; the effect of a filter can therefore be represented by an impulse 
response, which is the Fourier-transform of this gain factor, and will be denoted §(t) in the following. A filter used for 
the purpose of channel estimation is called a "Channel Estimation Filter" (CEF). The ability of a CEF to suppress noise 
has of course a direct effect on the performance of a coherent receiver. 

[0007] Among the parameters characterising a CEF, an important one is its bandwidth. Let us assume for the sake 

40 of the discussion that we have a "perfect" filter in the sense that its operation causes negligible distortion of the noiseless 
received signal. Clearly, this bandwidth should be large enough to avoid cutting a piece out of the signal waveform. 
On the other hand, if this bandwidth is too large, the filter will not eliminate as much noise as it should. 
[0008] A CEF ideally adapted to the statistics of a time-dispersive channel should be based on the so-called power 
delay spectrum of the CI R, which describes the expected channel output power for given delays. The major parameter 

45 of the power delay spectrum regarding the CEF is the above-mentioned excess delay of the channel, which represents 
the effective CIR length; namely, the expected power for delays greater than the excess delay is negligibly small. The 
problem is that this excess delay is unknown a priori to the receiver, and hence cannot be used in a straightforward 
manner to adapt the CEF bandwidth. Without adapting the filter, the CEF bandwidth has to be pre-set to a conservative 
value, according to the maximum expected delay in the planned applications of a communications system. This leads 

50 obviously to sub-optimal solutions, thus wasting some noise-suppression potential. 

[0009] Therefore, when the channel varies due to varying physical conditions, it is desirable to implement in the 
receiver a filter with adjustable filter parameters. Alternatively, if several filters for channel estimation are implemented 
in the receiver, and each is optimum for a different range of physical conditions, then a decision criterion is needed in 
order to select the one among these filters which will achieve the most efficient "de-noising" of the channel estimate. 

55 [0010] This problem arises, for example, in communications systems using Orthogonal Frequency-Division Multi- 
plexing (OFDM). 

[0011] OFDM is a multiple-channel modulation scheme. It is especially appropriate for highly frequency-selective 
channels such as typical channels for mobile communications, or for high-rate wireline transmission over copper lines. 
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Such channels are characterised by impulse responses which are substantially longer than one sample interval. This 
means that, even in the noiseless case, each received sample in the digital baseband domain is a superposition of 
multiple transmitted samples weighted by the appropriate channel coefficients. In order to resolve such "intersample 
interference", some kind of equalisation needs to be performed. 

5 [0012] The way OFDM combats intersample interference is by dividing the total channel bandwidth into a number 
D of substantially smaller portions, called subchannels. One OFDM channel comprises one parallel use of all subchan- 
nels. The data to be transmitted is collected into so-called OFDM symbols, and each OFDM symbol is transmitted in 
parallel on a number D u (0 < D u < D) of these subchannels. The transmitted subchannel signals are orthogonal to each 
other. Since the duration of one OFDM symbol is much longer than the sample interval, the problem of intersymbol 

10 interference is strongly reduced. 

[001 3] To get totally rid of intersymbol interference, a guard interval between two symbols is usually introduced during 
OFDM transmission. If the length of the guard interval exceeds the length of the channel impulse response, there is 
no residual intersymbol interference. Furthermore, if the guard interval is used in the form of a cyclic prefix, as is usually 
the case for OFDM transmission, one may implement a very simple equalisation of the frequency-selective channel in 

15 the frequency domain. 

[0014] The wireless Local Area Networks (WLAN) systems are examples of radio communications systems which 
use OFDM. ETSI (European Telecommunication Standard) BRAN (Broadband Radio Access Network) has developed 
a short-range high-data-rate communications system called "HIPERLAN type 2" (HI PERL AN/2). HIPERLAN/2 may be 
used to transport Internet Protocol (IP) packets, and will also be capable to act as a wireless Asynchronous Transfer 
20 Mode (ATM) system, as well as a public access system, e.g. with an interface to the Universal Mobile Telecommuni- 
cations System (UMTS). The physical layer of HIPERLAN/2 is based on OFDM, with a guard interval in the form of a 
cyclic prefix. Other WLAN systems based on OFDM have been standardised by ARIB in Japan (MMAC and its exten- 
sions), and IEEE in the US (IEEE802.11a and its extensions). 

[0015] In pilot-assisted channel estimation schemes for OFDM, known symbols are transmitted on given subchannels 
25 and time instants for training purposes; for example, in systems according to the HIPERLAN/2 or to the IEEE802.11a 
standard, there are two full OFDM pilot symbols preceding every burst of data-carrying OFDM symbols. The principle 
of decision-directed approaches for OFDM is quite similar: before channel estimation, some data symbols are being 
decided; these decided symbols are then treated in the same way as pilot symbols. 

[0016] The simplest channel estimation method consists in comparing a transmitted pilot symbol with the received 
30 value in the respective subchannel. The ratio between these two quantities then yields an estimated subchannel transfer 
factor. This method is known as the "least-squares estimation". 

[0017] One may alternatively estimate the channel while using a decision-directed coherent demodulation, as dis- 
closed in the paper by V. Mignone and A. Morello titled "CD3-OFDM: A Novel Demodulation Scheme for Fixed and 
Mobile Receivers", IEEE Trans. Comm., vol. 44, No. 9, pp. 1144-1151, Sep. 1996. In this method, received signals are 

35 fed back after they have been corrected by means of a channel coding scheme (which needs to be fairly powerful in 
order to ensure stability). Two types of appropriate filters are discussed: the first one has a pre-set bandwidth which 
is flat in the time-domain (and is of the order of the guard interval): such a filter is therefore not matched to the actual 
channel excess delay as it optimally should. The second one has a bandwidth which is adaptively determined after 
the receiver has calculated the CIR by Fourier-transforming the CTF: but this conversion from the frequency-domain 

40 to the time-domain is computationally expensive. 

[0018] A linear means for noise suppression in the described channel estimation process is a channel estimation 
filter exploiting correlations inside the channel transfer function, either in the frequency or in the time domain, or both. 
Assuming, for simplicity, that only correlations in the frequency domain shall be exploited in the OFDM signal, the 
characteristics of the channel transfer function over frequency can be regarded as a "band-limited" process. Band- 

45 limited means here that the transformation of the channel transfer function into the time-domain yields a CIR of limited 
length, i.e. the CIR length is substantially smaller than the OFDM symbol length. This usually holds for OFDM trans- 
mission systems. 

[001 9] Such an approach can be found in the paper by T. Onizawa et al. titled "A Simple Adaptive Channel Estimation 
Scheme for OFDM Systems", belonging to the same Proceedings of the 50 th IEEE Vehicular Technology Conference 

50 (VTC 1999-Fall), pp. 279-283, Amsterdam, The Netherlands, September 1999. In this paper, a local channel variation 
metric is used to select the appropriate CEF for a given portion of the entire bandwidth. For this purpose, the length 
of the difference vector between adjacent subchannels in the preamble is exploited to classify the frequency-local 
channel variation into three classes, each class being associated with a respective CEF. This method for exploiting 
the magnitude of difference signals, however, is derived only in a heuristic fashion, so that the level of performance of 

55 this approach is not reliable. 

[0020] Another example of a situation where adjusting the bandwidth of channel estimation filters is desirable con- 
cerns cellular systems based on Code-Division Multiple Access (CDMA). Wideband CDMA is used in the FDD (Fre- 
quency Division Duplex) mode of IMT 2000 (International Mobile Telecommunication)/UMTS (Universal Mobile Tele- 
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communication System). Also the International Standard I S-95, as well as evolutions thereof, are based on CDMA. 
[0021] A CDMA system may be viewed as an overlay of several spread spectrum (SS) signals. A SS modulator 
combines the data symbols with a "spreading sequence" (also called "signature sequence" or "spreading code"), at a 
sampling frequency (or "chip frequency") higher than the symbol frequency, and feeds the resulting sequence of chips 
5 into a pulse-shaping filter. The baseband signal generated this way is mixed up to the desired frequency band for 
transmission. 

[0022] In the receiver, the radio-frequency signal is mixed down to baseband and fed into a low-pass filter, which 
may be matched to the pulse-shaped chips. From the point of view of the receiver, the elements of the system situated 
before the pulse-shaping filter and after the chip-matched filter, taken together, can be considered as an "equivalent 
10 channel". 

[0023] A well-known receiver structure for SS is the "rake demodulator". A rake demodulator has several branches. 
Each respective branch is fed with the signal delayed by a so-called "tap delay", with the difference between the tap 
delays for any two successive branches being fixed. The signal samples, called "channel taps", are next de-spread by 
a so-called "correlator" in each branch. Finally, the outputs of all branches are combined according to well-known 
15 techniques. 

[0024] For more details about rake receivers, one may referto the book "Digital Communications" by John G. Proakis, 
Fourth Edition, McGraw-Hill International Edition, New York and Singapore, 2001. 

[0025] In the case of coherent demodulation, one usually applies so-called Maximum-Ratio Combining (MRC). This 
technique includes estimating channel taps of a reference signal, and using these estimates when combining all the 
20 branch outputs to obtain the demodulator output signal. 

[0026] In particular, in SS systems that use at least one pilot signal for channel estimation, the pilot signal has a 
spreading sequence different from the data signal. In each branch, there is, after tap delaying, one correlator for the 
data and one correlator for the pilot signal. 

[0027] A coarse estimate of the complex channel tap in a given branch is obtained at the output of the pilot correlator. 

25 However, this estimate usually includes a sizeable noise component. But since this noise has a large bandwidth, it is 
easily mitigated by a CEF. A CEF is therefore commonly provided for channel estimation in SS/CDMA receivers. If 
pilot symbols are transmitted in time-mutiplex with data symbols, one may associate with this CEF a functionality of 
interpolation with neighbouring data symbols. Moreover, a CEF reduces the impact of symbol decision errors in case 
one uses decision-directed channel estimation. 

30 [0028] Now, when the transmitter is on board of a vehicle, each channel tap may vary more or less randomly. This 
variation presents some patterns which can be modelled to some extent. It is caused by the movement of either the 
transmitter or the receiver, or also by the environment reflecting the transmitted signal. The evolution of each channel 
tap overtime can be regarded as a band-limited random process. It is important to note that the individual random 
processes for all branches might be different in general. The band limitation of each process is determined by the so- 

35 called Doppler bandwidth, or in general by the Doppler spectrum. It mainly depends on three terms, namely 

a) the velocity of the transmitter, the receiver, or the scattering environment, 

b) the RF channel frequency of the system, and 

c) the directions of branch signals emerging from the transmitter and impinging on the receiver. 

40 

[0029] A CEF ideally adapted to the execution-time channel characteristics should be based on this Doppler spec- 
trum. Since the Doppler spectrum is in general differentfrom one branch to another, an optimal CEF should be adapted, 
individually for every branch, to the actual Doppler characteristics of the respective branch. 

[0030] A major parameter characterising the Doppler spectrum is its bandwidth. Outside this bandwidth, the power 
45 contributions can be neglected. The optimum size of the CEF bandwidth should then be equal to the Doppler bandwidth 
of the considered branch. The problem here is that this Doppler bandwidth is unknown a priori to the receiver, and 
hence cannot be used in a straightforward manner to adapt the CEF bandwidth. 

[0031] In known rake receivers, the filters are usually simple low-pass filters. Thus, the CEF bandwidth has to be 
pre-set to a conservative value according to the maximum Doppler bandwidth, or, equivalently, to the maximum ex- 
50 pected velocity in all different possible application scenarios. But the range of mobile speed may be fairly wide, typically 
between 0 and 250 km/h, sometimes up to 500 km/h. This approach leads therefore in many cases to sub-optimal 
solutions where some noise-suppression potential is wasted. 

[0032] Generally speaking, one object of the present invention is to estimate the bandwidth of the noiseless compo- 
nent of any chosen quantity derived from a received signal (possibly the received signal itself), with the help of a set 
55 of alternative filters of known bandwidths (or, equivalently, with the help of a unique, adjustable, filter). 

[0033] Another object is to provide a method for selecting a filter out of a set of alternative filters (or, equivalently, 
for selecting the parameters of a unique, adjustable, filter) so as to obtain good noise-suppression from such a chosen 
quantity. 
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[0034] It is being assumed here that the bandwidth of the noise component of this chosen quantity is wider than the 
bandwidth of the noiseless component. 

[0035] Furthermore, the invention aims at fulfilling these objects by relying on "intrinsic" methods, viz., methods 
based on applying "in real time" these filters to the chosen quantity, and on analysing the output of the filters. 
5 [0036] Thus, according to a first aspect, the invention concerns, firstly, a method for bandwidth estimation in a data 
receiver, wherein said data receiver has been provided with a set of filters each having a known bandwidth and capable 
of removing noise from a quantity Z r derived from a received signal, the method comprising, for each one of these filters: 

a) the application of the filter to Z r , thus yielding a quantity Z F , 
10 b) the calculation of the power A of the difference between Z r and Z F , and 

c) the determination of whether the filter has truncated the noiseless component of Z r , on the basis of a criterion 
involving A and a 2 , where a 2 is an estimate of the power of the noise component of Z 

n n 

[0037] According to the same first aspect, the invention concerns, secondly, a method for adaptive filtering in a data 
15 receiver, wherein said data receiver has been provided with a set of filters capable of removing noise from a quantity 
Z r derived from a received signal, the method comprising: 

a) for each one of these filters: 

20 - the application of the filter to Z r , thus yielding a quantity Z F , and 

the calculation of the power A of the difference between Z r and Z F , and 

b) the selection of one filter out of the set on the basis of a criterion involving A and a 2 , where a 2 is an estimate 
of the power of the noise component of Z r , in such a way that the largest possible amount of noise is removed 

25 from Z r without truncating its noiseless component. 

[0038] The criterion according to the invention is based on the expected noise power which should be removed from 
the "noisy" quantity Z r for a given filter. The noise power in the channel is assumed to have been estimated independ- 
ently. 

30 [0039] The set of filters may contain either very many filters allowing for a very accurate bandwidth adaptation, or, 
for the sake of simplicity, only a few filters. Naturally, the invention may equivalently be implemented by adjusting the 
coefficients of a tuneable filter by means of successive trials. 

[0040] Compared to "extrinsic" adaptation methods which use a priori values (or, alternatively, require prior meas- 
urement operations), the present invention provides results derived from evaluations performed during the real-time 
35 functioning of the communications system. Furthermore, the method according to the invention requires a modest 
amount of computation, so that it is both convenient and accurate. 

[0041] In a preferred embodiment, said criterion requires that the power A be of the same order of magnitude as 
(a-a 2 ), where a is a number depending on the bandwidth of the respective filter and on the type of signal modulation 
boin n g used. 

40 [0042] Indeed, as shown below, this criterion leads to an optimal choice of filter in the case of several practically 
important types of signal modulations. 

[0043] The present invention may for example advantageously be applied to channel estimation in the case of a 
time-dispersive channel. 

[0044] In one such application, the type of signal modulation is OFDM, said quantity z r is a coarse estimate H coarse 
45 of the channel transfer factors vector, said quantity Z F is the filtered estimate H F of the channel transfer factors vector, 
and 



50 



where D is the total number of subchannels, D u the number of subchannels actually used, and D F the length of the 
bandwidth of the respective filter as discretised in the time domain. 

[0045] The invention allows in this case to estimate the excess delay of the CIR, and to achieve very good adaptive 
55 channel estimation by selecting the best CEF in the set in order to filter the CTF. 

[0046] It will be noted that the invention exploits the (in coherent OFDM receivers mostly already available) coarse 
and filtered estimates of the channel transfer factor. Only a small additional computational load is thus required to 
implement the invention in the receiver. 
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[0047] In another example, the present invention may advantageously be applied to channel estimation in the case 
of a transmission over a channel subjected to Doppler fading. 

[0048] In one such application, the type of signal modulation is Spread Spectrum (for example CDMA), the data 
receiver is a rake demodulator, said quantity Z r is a reference signal in a given branch of this rake demodulator, said 
5 quantity Z F is a channel tap estimate, and 

a=1 - B, 

10 where B is the equivalent noise bandwidth of the respective filter. 

[0049] The equivalent noise bandwidth of a filter is associated with the behaviour of its impulse response §(t) over 
time t; a precise definition will be given below. 

[0050] The invention allows in this case to estimate the Doppler bandwidth per branch, and to achieve very good 
adaptive channel estimation by selecting, for each channel tap, the best CEF in the set in order to filter the outcome 
15 of the pilot correlator. As an alternative complexity-saving solution, it is also possible to adapt the filter to the largest 
of all branch Doppler bandwidths. 

[0051] Clearly, for implementing the methods according to the invention, one needs to know the amount of noise 
contained in the received signal. The noise power can be evaluated by any known method. However, as a complement 
to the invention, the present disclosure will propose two methods for measuring noise power in systems where the 
20 transmitted signal includes a pilot signal, each method being adequate for a specific environment. 
[0052] Firstly, in SS systems, said noise power estimate a 2 may be obtained by: 

n 

high-pass or band-pass filtering of the pilot correlator output signal by means of a normalised filter (G) having a 
lower cut-off frequency at least as high as the maximum Doppler frequency that may occur at maximum transmitter 
25 speed, 

calculating the squared absolute value Isl 2 of the resulting filtered signal s, and 
applying an averaging filter (I) to Isl 2 . 

[0053] Secondly, if the pilot signal contains repeated sets of at least D amplitudes, a noise power estimate a 2 may 
30 be obtained from: 

$ 2 = 8-{l Vl -y 2 l 2 }/2D, 

35 where said sets of pilot amplitudes are being collected in D-dimensional vectors y, and y 1 and y 2 are two different 
occurrences of the pilot set, with e denoting expectation value. 

[0054] These methods for measuring the power of noise in the received signal are both reliable and easy to imple- 
ment. 

[0055] According to a second aspect, the invention concerns various devices. 
40 [0056] It thus concerns, firstly, a time-dispersion estimation device including: 

a set of channel estimation filters, each capable of filtering a previously-obtained coarse estimate of the channel 
transfer factors vector, 

a CTF-summation unit capable of calculating the difference between the coarse estimate and the filtered estimate 
45 of the channel transfer factors vectors, 

a stopband-power unit capable of calculating the expectation value A of the modulus-squared of said difference, and 
a criterion unit capable of calculating the value of 

so c = A-fc.a-a n , 

where a is a number depending on the bandwidth of the respective filter, a 2 is an estimate of the noise power in the 

received signal, and b is a predetermined number of the order of 1 . 

[0057] The invention concerns, secondly, a Doppler estimation device including: 

55 

a set of channel estimation filters, each capable of filtering a received pilot correlator output to yield a channel tap 
estimate, 

a tap-summation unit capable of calculating the difference between said pilot correlator output and said channel 
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tap estimate, 

a stopband-power unit capable of calculating the expectation value A m of the modulus-squared of said difference, 
and 

a criterion unit capable of calculating the value of 

5 

~ , A2 

C = A - £>-a-o n , 

where a is a number depending on the bandwidth of the respective filter, a 2 is an estimate of the noise power in the 
10 received signal, and b is a predetermined number of the order of 1 . 

[0058] The invention concerns, thirdly, a Doppler estimation device for a multiple-branch receiver, including, in each 
branch m (1 < m < L) : 

an identical set of channel estimation filters, wherein each filter is capable of filtering a received pilot correlator 
15 output to yield a channel tap estimate, 

a tap-summation unit capable of calculating the difference between said pilot correlator output and said channel 
tap estimate, and 

a stopband-power unit capable of calculating the expectation value A m of the modulus-squared of said difference, 
and including a criterion unit capable of calculating, when an identical filter is put in each branch m of the receiver, 
20 the value of 




where a is a number depending on the bandwidth of the respective filter, a is an estimate of the noise power in the 
received signal, and b is a predetermined number of the order of 1 . 
30 [0059] The invention concerns, fourthly, a noise-power estimation device, wherein the transmitted signal includes a 
pilot signal containing repeated sets of at least D amplitudes, including: 

a pilot-summation unit capable of calculating the difference between two received D-dimensional pilot vectors y 1 
and y 2 , and 

35 - a noise variance unit capable of calculating a noise power estimate according to: 

6n = e{l yi -y 2 l 2 }/2D, 

40 where 8 denotes expectation value. 
[0060] The invention also relates to: 

a modulated-signal reception apparatus including at least one device as succinctly described above, 
a telecommunications network including at least one reception apparatus as succinctly described above, 
45 - a data storage means containing computer program code instructions for executing steps of anyone method as 
succinctly described above, and 

a computer program containing instructions such that, when said program controls a programmable data process- 
ing device, said instructions mean that said data processing device implements anyone method as succinctly 
described above. 

50 

[0061] The advantages offered by these devices, apparatus, telecommunications networks, datastorage means and 
computer programs are essentially the same as those offered by the methods according to the invention. 
[0062] Other aspects and advantages of the invention will emerge from a reading of the detailed description, which 
will be found below, of preferred embodiments given by way of non-limitative example. This description refers to the 
55 accompanying drawings, in which: 

Figure 1 is a diagram illustrating the standard time-domain structure of an OFDM preamble, 

Figure 2 depicts schematically a time-dispersion estimation device according to a first embodiment of the invention, 
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and also a suitable noise power estimator, 

Figure 3 depicts schematically the conventional functioning of a rake demodulator, 

Figure 4 depicts schematically a Doppler estimation device according to a second embodiment of the invention, 
and also a suitable noise power estimator, and 
5 - Figure 5 is a flow diagram illustrating the main steps leading to the adaptive selection of an optimal channel esti- 
mation filter according to either the first or the second embodiment of the invention. 

[0063] Two preferred embodiments of the invention will be described in detail below. The first embodiment concerns 
OFDM modulation schemes, and the second embodiment concerns SS/CDMA communications systems. 

10 

FIRST EMBODIMENT OF THE INVENTION 

[0064] For introduction, consider the transmission of one single OFDM symbol generated with a D-dimensional In- 
verse Discrete Fourier Transform (IDFT). Each of the D subchannels is modulated by a complex-valued amplitude 
15 [u], where v is the subchannel number (0 < v < D) representing subchannel frequency, and \i indicates the OFDM 
symbol number representing transmission time. The data are mapped on the subchannel amplitudes A^[v]. The OFDM 
symbol is equipped with a guard interval of D g samples. The discrete-time complex baseband OFDM transmitted signal 
for OFDM symbol number \i, when written out as a vector over discretised timep, with: - D g <p < D, has the following 
components: 

20 
25 

which can be generated by an IDFT of dimension D. 

[0065] Using e to denote expectation value (computed over a set of samples), the average signal power is 

-MkH| 2 }=f ■ 

where E s is the average energy per channel symbol and Tis the modulation interval. 

35 [0066] For convenient notation, the symbol index [i, will be suppressed in the following. 

[0067] We assume transmission over a multipath channel, whose CIR will be modelled, during the transmission of 
each OFDM symbol, by means of a discrete-time vector h having components h[p] (p = 0,..., D-1). The fact that the 
continuous-time CIR is non-zero only up to the excess delay t is equivalent to having only the components h[p] with 
p = 0,..., D e -1 (where D e is the discretised excess delay of the CIR) being non-zero. 

40 [0068] As discussed in the introduction, the noiseless received signal is given by 

y=a*h, 

45 where * denotes circular convolution. 

[0069] For convenient notation, without loss of generality, perfect frequency synchronisation will be assumed in the 
receiver. At the receiver input, samples n[p] of Additive White Gaussian Noise (AWGN) corrupt the noiseless signal 
samples y [p] to yield the actually-received samples 

50 

y[p] = y[p] + n[pi 

[0070] The noise power is 

55 

= e {ln[p]l 2 } = ^ , 
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where N 0 is the "one-sided power spectral density" of the white noise. The channel Signal-to-Noise power Ratio (SNR) 
at the receiver input is thus 

E s /N 0 = o 2 5 /o 2 n . 

[0071] Implementing a Discrete Fourier Transform (DFT), one defines the discrete Channel Transfer Factor (CTF) 
at subchannel v as 



D-1 



H\o]m^h\p]e ° 

p=0 



2% 



(2) 



[0072] The OFDM receiver yields the noisy subchannel amplitudes via DFT processing as 



^ D-1 

5 p. 

= Y[v] + N[v] , 



2k 



where 



y[u] = 4u] H l»] 



is the noiseless received subchannel amplitude, and N[v] is the resulting noise sample on subchannel v. Here we 
assumed that the guard interval D g is sufficiently large so that the linear convolution of the transmitted OFDM signal 
and the CIR is equal to their circular convolution in the time interval of interest. 

[0073] It is obvious from Equation (3) that for data recovery from Y[v], an estimate of /-/[d] is required. 
[0074] We shall now describe, as an example of a method which can be implemented in the receiver to obtain a 
coarse estimate of H[d] within the framework of the present invention, the known "least-squares" estimation method. 
[0075] The active subchannel positions will be referred to by the subchannel indices v which belong to a set U 
containing the D u indices of these positions. The remaining (D - D u ) subchannels, namely the ones in which the mod- 
ulation amplitude is set to zero, are usually called "virtual" subchannels. 

[0076] We shall assume that, for so-called training purposes, a preamble symbol containing known subchannel am- 
plitudes covering the entire active frequency-domain is provided; we will collect these amplitudes in a vector A p of 
length D u . 

[0077] ) The transmitted time-domain preamble (or pilot) OFDM symbol is depicted in the upper part of Figure 1 . 
The repeatedly transmitted time-domain sequence of length D is obtained by the I DFT 



(4) 



where we introduced the D u xD partial DFT matrix 



2te 
- J — up 
D 



(5) 



pejp,...,D-l} 
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and ^denotes the conjugate transpose of a matrix. 

[0078] The sample sequence resulting from the convolution of the transmitted preamble sample sequence with the 
CIR h is disturbed by additive noise. Since the guard interval and the repeated preamble, which are assumed to be of 
sufficient size, are both cyclic, the linear convolution is equal to the circular convolution for these signal segments 
which are of interest for final demodulation. Thus, the noiseless received sequence 9 p , of length D is equal to the 
circular convolution of a p and h. 

[0079] In the lower part of Figure 1 , the temporal structure of the received signal is indicated. Perfect channel fre- 
quency synchronisation has been assumed. The successive received sequences y-, y 2 , and so on, are made of the 
repeated noiseless sequence y p disturbed respectively by the (statistically independent) successive D-dimensional 
noise vectors n 1 , n 2 , and so on. 
[0080] Implementing a DFT, we find: 

Y p = -LD u y p =diag{-L D u a p }• D u h = diag { A p } H = A H , (6) 

VZ> vD 

where we collected the D u channel transfer factors h^\>] (v GU) in the column vector 

H = D u h , (7) 
and we collected the transmitted pilot subchannel amplitudes in the D u x D u diagonal matrix 

A - diag{A p }. (8) 

[0081] Each received sequence is of the form: 

y P = y P + n , (9) 

where n is a noise vector. 

[0082] This received signal y p is transformed into frequency domain via the DFT 



Y p - j=D u y p = A-H + -j=D u n = A-H + N . (1 0) 



where we used Equation (6), and we introduced the D u -dimensional frequency-domain noise vector N. 

[0083] Finally, the "least-squares estimate" H LS of the CTF vector is obtained by dividing out the frequency-domain 

received signal by the pilot subchannel amplitudes: 

H LS =A 1 Y P . (11) 



[0084] This operation can be seen as extracting or removing the influence of the transmitted pilot amplitudes from 
the received signal. 

[0085] The least-squares estimation method just described is only one possible way to obtain a coarse estimate of 
the CTF. Once a coarse estimate H coarse has been obtained, filtering of H coarse by the CEF yields the filtered estimate 
H F of the CTF vector. The filters used in the framework of the present invention are assumed to be normalised in such 
a way thatthere is neither an amplification nor an attenuation in the passband; in otherwords, onetries, when designing 
a filter, to have the parts in the filter input signal which belong to the filter passband be equal to the corresponding parts 
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in the filter output signal, up to signal distortions introduced by the filter itself. Then, a "good" filter is a filter which, 
besides introducing negligible distortions, is characterised by a passband which is well-adjusted to the noiseless signal 
characteristics. 

[0086] We shall now describe a method according to the invention for estimating the effective length D e of the CIR 
on the basis of the measurement of power in the stopband subspace of the^filter. Shortly speaking, one obtains to 
begin with, for example in the above-described manner, a coarse estimate H coarse of the CTF vector: this estimate 
plays here the role of a quantity Z r representing the received signal. Several different channel estimation filters, each 
having a discretised bandwidth of Respective length D F , are then applied to Z r in order to obtain respective quantities 
Z F , namely the filtered estimates H^. Finally, one computes^he expectation value (averaged over frequency) of the 
difference between the filter "input" H coarse , and the "output" H F of each filter. The derivation below leads to a decision 
criterion for choosing the best filter in the set on the basis of the values thus computed, and shows by the same token 
how an estimate for the CIR length can be obtained. 

[0087] Consider, for any given filter, the expectation value A of the difference between the "coarse" and the "filtered" 
estimates of the CTF: 



H coarse ^1 



(12) 



A 

where averaging overthe D u components of the CTF vectors is implied. Since H F results, by definition, from the filtering 
of H coarse this quantity A can be interpreted as the average power in the stopband subspace of the CEF. 
[0088] Furthermore, A can be written as: 

A=P n +P t .(13) (13) 

Indeed, the power A consists of two components. The first component is the noise power P n in the stopband subspace 
of the filter, which is thus removed thanks to the filter. The second component is the tail power P t in the portion of the 
CIR which has been truncated by the CEF in the case where the length D F of the filter is shorter than the excess delay 
D e of the CIR; in order to make P t zero, as it should optimally, one needs then to make D F larger than D e . However, 
making D F too large results in an inefficient noise removal. 

[0089] We shall assume that these two power components are uncorrelated, and also that the components h[p] of 
the CIR are uncorrelated. We shall also assume that the CEF is "perfect" in the sense that it is able to extract the D F 
components of the CIR which are located in the bandwidth of the CEF without any distortion. Also, the noise afflicting 
the received signal will be assumed to be white noise. Finally, the pilot subchannel amplitudes will be assumed to 
correspond to points on the unit circle, i.e., AA+ = I. 
[0090] After some calculus, one finds then: 

Pn»2-°n%-(D-D F ), (14) 



and 

p=£)p+1 

[0091] Therefore, taking, for the purpose of this first embodiment: 



a = ^(D-D F ),(16) (16) 



the sought-for criterion reads: 
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• if: A - a-a 2 > 0 , 
then: D~ > D F , and 

• if: A-x-o =0, 
then: D e <D F . 

5 

[0092] In practice, the value of the noise power must be measured in one way or another, and the criterion as im- 
plemented will actually use the value a 2 resulting from this measurement. 

[0093] The above criterion has been derived on the basis of a few idealistic assumptions. Now taking practical im- 
pairments into account, the time-dispersion classification should preferably be defined as follows: 

10 

• if: C> q, then: D e > D F , and (17a) 

• if: C< q, then: D e < D F , (17b) 

where: 

15 

C-A-fc-a-a 2 ,. (17c) 

In Equations (17a-c), the factor b accounts for the filter not being perfect, and q can be used to influence the threshold 
20 decision. Usually, 1 < b< 2 and q~ 0 are reasonable values. The factor b is required due to the perfect-filter assumption, 
which disregards the passband distortion by the filter, as well as the imperfect stopband behaviour, and the non-ideal 
transition region between pass- and stopband. Furthermore, the filter transients in the proximity of subchannels without 
an available pilot are neglected in theory. The imperfect filter characteristic leads to an actual value for A which is higher 
than is expected in our theory, so that the theoretical difference power needs to be increased by this factor b. 
25 [0094] An exemplary device (1 00) for implementing this first embodiment of the invention is schematically depicted 
in Figure 2. 

[0095] A received pilot sequence is Fourier-transformed, according to Equation (1 0), in unit (1 02), after which a least- 
squares estimate of the CTF vector is obtained by "removing pilot amplitudes", accordingto Equation (11 ), in unit(103). 
[0096] The following steps are implemented for each filter in a set (1 04) of channel estimation filters. First, the least- 

30 squares estimate of the CTF vector is filtered by this filter. The difference between these two estimates is calculated 
in CTF-summation unit (105), after which the power A of this difference is calculated in stopband-power unit (106). 
Finally, after reading the value of a for this particular filter and an estimate for the noise power in the received signal 
provided by noise-power estimator (200), criterion unit (107) calculates the value of C associated with this filter. 
[0097] A method, based on this criterion, leading to the selection of the optimal CEF will be described below in 

35 reference to Figure 5. 

[0098] A method to estimate said noise power in the present context will now be described. 

[0099] Let us refer again to Figure 1 . Exploiting the fact that the useful signal portion is identical in the two received 
signal parts of the repeated preamble, and using previous assumptions and definitions, one derives easily: 

8{l yi -y 2 U2D$ 2 . 
[0100] Hence, a noise variance estimate is obtained by 

6n=e{l yi -y 2 l 2 }/2D,(18) (18) 
where the expectation value is computed in the time domain. 

[0101] Clearly, it is not compulsory to use the entire vectors y., and y 2 ; any subvectors of y-, and y 2 could actually 

50 be used in order to estimate the noise variance. 

[0102] It is to be noted that this estimator is not limited to OFDM preambles; it can as well be used with repeated 
single-channel preambles. Generally speaking, applications of this noise power estimator are not limited to its combi- 
nation with the channel estimation according to the invention, the only requirement for its applicability being the exist- 
ence of a repetition preamble. For example, it can also be used for link-quality control and adaptation. 

55 [0103] Going back to Figure 2, one finds, correspondingly, a schematic illustration of a device (200) for noise power 
estimation which can be implemented in a communications system using a repetition preamble. The difference between 
two received preamble vectors is calculated in pilot-summation unit (201), and is used in noise variance unit (202) to 
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calculate an estimate for the noise power in the received signal, according to Equation (18). In the present context, 
this noise estimate is finally passed on to criterion unit (107) for the purpose of calculating the criterion according to 
the invention, as described above. 

[0104] It should be noted that the existence of a repetition preamble is here required only when implementing this 
5 noise power estimator. If a noise variance estimate is obtainable from some other source, then the time-dispersion 
estimator according to the invention does not require a repetition preamble. 

SECOND EMBODIMENT OF THE INVENTION 

10 [0105] Figure 3 depicts schematically the conventional functioning of a rake demodulator, in the case of SS systems 
that use at least one pilot signal for channel estimation. 

[0106] After being filtered by a chip-matched filter in unit (301), the baseband received signal is dispatched between 
the different branches of the rake receiver. The signal in each branch m is suitably delayed in unit (302). Then data 
and pilot symbols are sent to their respective correlators (303) and (304). 
15 [0107] The output of the pilot correlator (304) is made of symbols separated by time-intervals T p . Each of these 
symbols can be written as 



Xm=9m+"m>V*) ( 19 ) 

where g m and n m denote the ideally-noiseless channel tap, and the noise (assumed to be AWGN) respectively. For 
convenience, the time index is suppressed here and in the following. 

[0108] The filtering of x m by a low-pass filter in unit (305) yields the channel tap estimate g m , which can be expressed 
as the convolution of x m with the impulse response $ of this CEF: 

9m=9m**+"m**- ( 20 ) ( 20 ) 



[0109] The complex conjugate of this estimate is calculated in unit (306) and is multiplied, in unit (308), by the output 
of the data correlator (303) after the latter has preferably received compensation, in unit (307), of the group delay 
caused by the CEF. Finally, the resulting signals in all branches are summed up in unit (309). 

[0110] We shall now present a derivation leading to the criterion for filter selection according to this second embod- 
iment of the invention. 

[01 1 1 ] This selection is performed after symbol x m , which plays here the role of a quantity Z r representing the received 
signal, has come out of the pilot correlator (304), and before the corresponding channel tap estimate g m , which plays 
here the role of a quantity Z F resulting from the filtering of Z r , enters unit (306). 

[01 1 2] Consider then the expectation value A m of the difference between the input x m and the output g m of the channel 
estimation filter: 



40 



A, 



m 



■Si 



m\ 



(21) 



As in the first embodiment of the invention, this quantity A can be interpreted as the average power in the stopband 
subspace of the CEF. 

[0113] Using Equations (19) and (20), one gets 



50 



*m=s\\g m -gm*$ 2 \ + *\\ n m- n m*$\ 2 \ ■ < 22 ) 



[0114] With the assumption of the noise being white noise, the last term in Equation (22) is simply calculated as 

55 
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5 



n -n 
m m 



4>| 2 j = c^(l-i?) , (23) 



where a 2 denotes the power of the noise n m , and B is the "equivalent noise bandwidth" of the filter, which can be 
defined as 

10 B = pW , (24) 

IZ p *[p]l 2 

and is always comprised between 0 and 1 . It is assumed here, as in the first embodiment, that the filter response is 
normalised, so that the input signal is neither attenuated nor amplified in the passband of the filter. 
15 [0115] The term e{lg m - g m * $l 2 } can be viewed as a measure of the distortion of the noiseless signal waveform by 
the filter. If we were dealing with a filter having constant magnitude response and constant group delay in its passband 
up to its cut-off frequency, this cut-off frequency would be equal to (B/2 T p ). In this ideal case: 

2 o e {\g m - g m * *l 2 } = 0 for (B/2T p ) > f D , (25a) 

and 

25 *V9m- 9m ^> 0 for ( B/2T p) < f D> ( 25b ) 

where f D denotes the maximum Doppler frequency (corresponding to a Doppler bandwidth of 2f D ). 

[0116] In practice, the best filter in the available set will therefore be the one having the smallest bandwidth among 

those who verify: 

30 

C<q, (26a) 

where 



35 



and 

40 



C = A m -b-a-S 2 n , (26b) 



a = 1 - B (26c) 

The parameters b and q appearing here are analogous to the ones introduced in the first embodiment. It should be 
45 noted that C may become negative due to estimation errors. 

[0117] An exemplary device (400) for implementing this second embodiment of the invention is schematically de- 
picted in Figure 4. 

[0118] The following steps are implemented for each filter in the set. First, the difference between X m and g m is 
calculated in tap-summation unit (401), after which the power A m of this difference is calculated in stopband-power 
50 unit (402). Finally, after reading the value of this power and of the parameter a for this particular filter, and an estimate 
for the noise power in the received signal provided by noise-power estimator (500), criterion unit (403) calculates the 
value of C associated with this filter. 

[0119] The filter selection method described so far evaluates the signals in any single branch of a receiver. It may 
however be reasonable, when the receiver, for example a rake receiver, has multiple branches, to assume the same 
55 Doppler bandwidth as well as noise spectrum in all branches. The accuracy of the filter selection scheme can, therefore, 
be increased by considering all branches in a parallel computation , assuming one puts an identical filter in each branch 
of the receiver. In this case, preferably, the respective value of the power A m for each of the L branches will be input 
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to criterion unit (403), and the decision variableto be evaluated in the Doppler estimation device (400') for each available 
type of filter will be: 



10 rather than Equation (26b). 

[0120] The present embodiment of the invention can be applied to any SS communications scheme, even if it does 
not use actual CDMA for user multiple access. This is, for instance, true for I EEE802. 1 1 , where SS is used with channel- 
sense multiple access, but without CDMA in the sense of simultaneous overlay of SS signals. This embodiment can 
also be applied to any SS system that uses channel estimation filters, even if it does not use training symbols; for 

15 instance, it may be used with decision-directed channel estimation. 

[0121] Finally, it should be noted that even though, for convenience sake, the above discussion dealt with a rake 
receiver without interference cancellation, the present invention is not restricted to such a receiver. Generally speaking, 
this embodiment of the invention exploits the temporal correlations in a channel estimate in order to achieve improved 
noise-suppression for channel estimation. Any kind of post-processing that makes use of a channel estimate can thus 

20 benefit from the invention. The standard linear rake receiver (having combining weights equal to the complex conjugate 
of the tap weights) was only presented above as one example of how the invention may be usefully implemented. 
[0122] A method, based on the above-derived criterion, leading to the selection of the optimal CEF will be described 
below in reference to Figure 5. This method for adaptive channel estimation is based on the use of a multiplicity of low- 
pass filters, whose bandwidths are predetermined and scaled by reasonable steps, so as to cover the whole range of 

25 possible Doppler bandwidths. Actually, considerable improvement can already be expected by using just two or three 
alternative filters; for example, for a system with maximum Doppler bandwidth of 300 Hz, a choice of filter bandwidths 
with three alternative filters might be: 300 Hz, 150 Hz and 50 Hz. 

[0123] This second embodiment of the invention requires a noise power estimate. Going back to Figure 4, one finds 
a schematic illustration of a device (500) for noise power estimation in a band-limited signal. 
30 [0124] First, the pilot correlator output signal x m is filtered by means of a high- or band-pass filter (G) in unit (501) (it 
will be assumed below that this filter is normalised so that the sum of the absolute value squared of the filter coefficients 
is equal to one). Filter (G) shall have a lower cut-off frequency at least as high as the maximum Doppler frequency that 
may occur at maximum mobile speed. 

[0125] The absolute value squared of the resulting signal is next calculated in unit (502), and submitted to an aver- 
ts aging filter (I) in unit (503). The noise power estimate o 2 for the symbol-spaced signal is then equal to the output of 
this averaging filter (I). 

[0126] The accuracy of this noise power estimate can obviously be improved, for a multiple-branch receiver, by 
averaging the branch noise power estimates over several, or all branches. 

40 FILTER SELECTION ALGORITHM 

[0127] Figure 5 is a flow diagram illustrating the main steps of an exemplary process leading to the selection of an 
optimal channel estimation filter according to the invention. For simplicity sake, we shall consider only the case (actually 
rather general) where the optimal CEF is such that the power A is of the order of (a-a 2 ), so that the criterion of Equation 
45 (26a), with the appropriate values of a, b and q, is applicable. 

[0128] The receiver has at its disposal a set of k filters, which will conveniently be labelled in the order of increasing 
filter bandwidth. For example, in the above-described embodiments of the invention, this means increasing value of 
D F for the first embodiment, and increasing values of B for the second one. 

[0129] The selection process is initialised at step 601 by assigning the value zero to a counter /. This counter is 
50 incremented at step 602 and will serve in the following as filter label. During step 603, the value C, of Cfor filter / is 
calculated using appropriate methods; for example in the case of the above-described embodiments, one may use the 
methods described in reference to Figure 2 and Figure 4 respectively. During step 604, this number C f is compared to 
a predetermined number q. If C,is larger than q, and if, as checked at step 605, there remains at least one yet untested 
filter in the available set, then the filter under consideration is considered as not optimal, and one goes back to step 
55 602. If, on the other hand, either C, is smaller than q (or equal to q within the precision of the computation), or the filter 
under consideration is the last one, then, during step 606, this filter is taken to be the optimal filter and the process is 
stopped. 
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A method for bandwidth estimation in a data receiver, wherein said data receiver has been provided with a set of 
filters each having a known bandwidth and capable of removing noise from a quantity Z r derived from a received 
signal, the method comprising, for each one of these filters: 

a) the application of the filter to Z r , thus yielding a quantity Z F , 

b) the calculation of the power A of the difference between Z r and Z F , and 

c) the determination of whether the filter truncated the noiseless component of Z r , on the basis of a criterion 
involving A and a 2 , where a 2 is an estimate of the power of the noise component of Z r . 

A method for adaptive filtering in a data receiver, wherein said data receiver has been provided with a set of filters 
capable of removing noise from a quantity Z r derived from a received signal, the method comprising: 

a) for each one of these filters: 

the application of the filter to Z r , thus yielding a quantity Z F , and 

the calculation of the power A of the difference between Z r and Z F , and 

b) the selection of one filter out of the set on the basis of a criterion involving A and a 2 , where a 2 is an estimate 
of the power of the noise component of z r , in such a way that the largest possible amount of noise is removed 
from z r without truncating its noiseless component. 

A method according to Claim 1 or Claim 2, characterised in that said criterion requires that the power A be of the 
same order of magnitude as (a-o 2 n ), where a is a number depending on the bandwidth of the respective filter and 
on the type of signal modulation being used. 

A method according to Claim 3, characterised in that said type of signal modulation is OFDM, said quantity Z r is 
a coarse estimate H coarse of the channel transfer factors vector, said quantity Z F is the filtered estimate H F of the 
channel transfer factors vector, and 

«=2B<0-z*)- 

where D is the total number of subchannels, D u the number of subchannels actually used, and D F the length of 
the bandwidth of the respective filter as discretised in the time domain. 

A method according to Claim 3, characterised in that said type of signal modulation is Spread Spectrum, the 
data receiver is a rake demodulator, said quantity Z r is a reference signal in a given branch of this rake demodulator, 
said quantity Z F is a channel tap estimate, and 

a=1-B, 

where B is the equivalent noise bandwidth of the respective filter. 

A method according to Claim 5, wherein the transmitted signal includes a pilot signal, characterised in that said 
noise power estimate a 2 is obtained by: 

n 

high-pass or band-pass filtering of the pilot correlator output signal by means of a normalised filter (G) having 
a lower cut-off frequency at least as high as the maximum Doppler frequency that may occur at maximum 
transmitter speed, 

calculating the squared absolute value Isl 2 of the resulting filtered signal s, and 
applying an averaging filter (I) to Isl 2 . 

A method for measuring noise power in a received signal, wherein the transmitted signal includes a pilot signal 
containing repeated sets of at least D amplitudes, characterised in that a noise power estimate a 2 is obtained 
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from: 

° 2 n = zVy^ -y 2 |2 } / 2D ) 

where said sets of pilot amplitudes are being collected in D-dimensional vectors y, and y-| and y 2 are two different 
occurrences of the pilot set, with 8 denoting expectation value. 

. A time-dispersion estimation device (100), including: 

a set (1 04) of channel estimation filters, each capable of filtering a previously-obtained coarse estimate of the 
channel transfer factors vector, 

a CTF-summation unit (1 05) capable of calculating the difference between the coarse estimate and the filtered 
estimate of the channel transfer factors vectors, 

a stopband-power unit (106) capable of calculating the expectation value A of the modulus-squared of said 
difference, and 

a criterion unit (1 07) capable of calculating the value of 

C = A-b-a-o*, 

where a is a number depending on the bandwidth of the respective filter, a 2 is an estimate of the noise power in 
the received signal, and b is a predetermined number of the order of 1 . n 

. A Doppler estimation device (400), including: 

a set (305) of channel estimation filters, each capable of filtering a received pilot correlator output to yield a 
channel tap estimate, 

a tap-summation unit (401 ) capable of calculating the difference between said pilot correlator output and said 
channel tap estimate, 

a stopband-power unit (402) capable of calculating the expectation value A of the modulus-squared of said 
difference, and 

a criterion unit (403) capable of calculating the value of 

C = A-b-a-o n , 

where a is a number depending on the bandwidth of the respective filter, a 2 is an estimate of the noise power in 
the received signal, and b is a predetermined number of the order of 1 . 

0. A Doppler estimation device (400 1 ) for a multiple-branch receiver, including, in each branch m (1< m < L): 

an identical set (305) of channel estimation filters, wherein each filter is capable of filtering a received pilot 
correlator output to yield a channel tap estimate, 

a tap-summation unit (401 ) capable of calculating the difference between said pilot correlator output and said 
channel tap estimate, and 

a stopband-power unit (402) capable of calculating the expectation value A m of the modulus-squared of said 
difference, and including a criterion unit (403) capable of calculating, when an identical filter is put in each 
branch m of the receiver, the value of 



1 2 

L m=\ 

where a is a number depending on the bandwidth of the respective filter, o 2 is an estimate of the noise power in 
the received signal, and b is a predetermined number of the order of 1 . 
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11 . A noise-power estimation device (200), wherein the transmitted signal includes a pilot signal containing repeated 
sets of at least D amplitudes, including: 

a pilot-summation unit (201) capable of calculating the difference between two received D-dimensional pilot 
vectors y 1 and y 2 , and 

a noise variance unit (202) capable of calculating a noise power estimate according to: 

6^8{l yi -y2l 2 }/2D, 

where e denotes expectation value. 

12. A modulated-signal reception apparatus, characterised in that it includes at least one device according to Claim 
8 and/or Claim 9 and/or Claim 1 0 and/or Claim 1 1 . 

13. A telecommunications network, characterised in that it includes at least one reception apparatus according to 
Claim 12. 

14. A data storage means, characterised in that it contains computer program code instructions for executing steps 
of a method according to any one of Claims 1 to 7. 

15. A data storage means according to Claim 14, characterised in that it is partially or totally removable. 

16. A computer program, characterised in that it contains instructions such that, when said program controls a pro- 
grammable data processing device, said instructions mean that said data processing device implements a method 
according to any one of Claims 1 to 7. 
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